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FINITE ELEMENT STRESS ANALYSIS OF 
IDEALIZED COMPOSITE DAMAGE ZONES 

ABSTRACT 

A quasi three dimensional finite element stress analysis of 
idealized damage zones in composite laminates is presented. The 
damage zones consist of a long centered groove or cutout extending 
one or two layers in depth from both top and bottom surfaces of a 
thin composite laminate. 

Elastic results are presented for compressive loading of four and 
eight layer laminates. It is shown that a boundary layer exists near 
the cutout edge similar to that previously shown to exist along free 
edges. The cutout is shown to produce significant interlaminar stresses 
in the interior of the laminate away from free and cutout edges. The 
interlaminar stresses are also shown to contribute to failure which is 
defined using the Tsai-Wu failure criteria. 

Nonlinear stress-strain curves are presented for compressive 
loading of three configurations of a [90/±45/0]^ laminate. It is 
shown that a damaged laminate fails at approximately 60 percent of 
load carried by an undamaged laminate. It is also shown that initial 
failure occurs at the free edge in the 0® layer. Subsequent failure 
in the interior of the 0® layer, as well as at the free and cutout 
edges of the other layers is shown to lead to total failure. 

The material system studied is a o**aphite-epoxy, Thornel 300/ 

Narmco 5208, with nonlinear material behavior represented by 
Ramberg-Osgood approximations. 
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Chapter 1 
INTRODUCTION 

The expanding use of laminated fibrous composi .es in such struc* 
tures as aircraft ailerons and turbine blades has prompted many studies 
into impact damage. The American Society for Testing and Materials has 
held a symposium on Foreign Object Damage (FOO) and a *^ent survey at 
the NASA/Langley Technical Library lists 272 separate articles on impact 
in composites. Thus, impact damage has became a major concern to the 
design engineer. 

The major emphasis, though, has been on failure mechanisms, ii^>act 
response and residual strength. Further, most studies have been con- 
cerned with damaged composites where the damage extends through the 
total thickness, such as holes. Little attention has been given to 
surface damage which "peels" off only a layer or two of the laminated 
conposite. Further, few studies have been made into the stress distri- 
bution around the damaged region or the effect of the free edge at the 
damaged region on the three-dimensional stress distribution. 

This study examines an idealized damaged region: a centered rec- 

tangular groove or cutout, extending one or two layers in dep end 
running the entire length of a long thin symmetric composite laminate. 
For synwerry reasons, the tout is on both the top and bottom surfaces 
(see Fig. 1). This idealized damaged region, which simulates several 
layers "peeled" off, was chosen as a tractable problem providing insight 
into the more general problem of the stress distribution near an ir- 
regularly shaped damage zone. 


1 



2 


The idealized damage zone is studied by the finite element method. 
The finite element programs NONCOM, developed by Renieri and 
Herakovich [1], and Herakovich, Renieri and Brinson [2], and a modified 
version, NONCQMl, developed by Humphreys and Herakovich [3], are both a 
quasi-three dimensional analysis. Hygrothermal and nonlinear capa- 
bilities are included in the programs, though this study will not 
examine hygrothermal effects. Some modifications have been made to 
NONCQNl to improve efficiency and running time and the modified version 
is called N0NC0M2. In addition, a Tsai-Hu failure criteria has been 
included and used to predict failure. 

Due to the increased concern over the compressive behavior of 
laminated composites, all loadings are conq>ressive. A number of lam- 
inates are studied, with and without cutouts, including cross-ply, 
angle-ply, and quasi-isotropic laminates. The material system examined 
is Thomel SOO/.iarmcc 5208 graphite-epoxy for which a large body of data 
exists in the literature. 



Chapter 2 
LITERATURE REVIEW 

Idealized damage zones have a long history of investigation. 
Timoshenko and Woinowsky-Krieger [4] developed analytic expressions for 
the stress distribution near a circular hole in a uniformly loaded 
infinite plate. Timoshenko and Goodier [5] extended the analysis to 
plates with elliptical holes and fillets. Both references used a plane 
strain, linear elastic, isotropic approach and neglected o^, and 

presented stress concentration values for the 

cases examined. 

Roark and Young [6], as well as Peterson [7], give stress concentra- 
tion factors for various material discontinuities such as holes, grooves 
and notches. The stress concentration factors are merely con4>11ed in 
references [6] and [7], having been experimentally determined by other 
authors. The data presented were for linear elastic. Isotropic 
materials only. None of the discontinuities presented had the same 
configuration as the cutout examined in thU study. 

Theoretical stress distributions around holes In anisotropic plates 
have been given by Savin [8]. Circular, elliptical and filleted 
rectangular holes are studied using a complex mapping approach. Again, 
a plane strain, linear elastic approach Is utilized and o^, and 
stresses are neglected. 

Experimental studies into the stress state near cutouts in composite 
plates have been conducted by Rowlands, Daniel and Whiteside [9]. Cir- 
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cular, elliptical and square cutouts in boron epoxy composite plates 
were investigated. Two laminate configurations and a number of hole 
diameter-to-width ratios and hole diameter-to-thickness ratios were 
studied. Another s :dy by Ashton, Burdorf and Olson [10] gave experi- 
mentally determined stress concentration values for cutouts in graphite- 
epoxy coupons. Both of these studies were concerned with the hole's 
effect on the strength of the laminate and neither investigated inter- 
laminar stresses around the hole. 

A three dimensional finite element program developed by Barker, 

Dana and Pryor [11] was used to analyze cutouts in finite width laminates 
subjected to uniform axial strain. Interlaminar stresses (o , x and 

Z yZ 

cutouts were investigated and stress concentration values 
for circular, square and diamond-shaped holes were gi» Nonlinear 
effects were not considered, nor was failure of the laminate studied. 

The idealized damage zone examined in this study has apparently not 
been investigated before. A major reason for this nay be the lack of 
analytic tools to handle this problem. The development of the finite 
element program NONCOM in references [1,2] has provided the necessary 
analytical tools for such an analysis- 



Chapter 3 

THEORETICAL BACKGROUND 

The problem under consideration is the stress analysis of a long, 
finite width, synmetric composite laminate with a centered strip removed 
along the entire length of the laminate (Fig. 1). The laminates may be 
subjected to uniform thermal, moisture, and axial mechanical loading and 
all stresses and strains are assumed to be independent of the axial (x) 
coordinate. The resulting differential equations were presented by Hsu 
and Herakovich [12] while the theory behind the finite element program 
N0NC0N2, tdtich obtains a solution to the problem defined by those differ- 
ential equations, was presented in *%ferences [1*2]. This report will 
present only the main points of the theoretical development. 

3.1 Problem Formulation 

The laminate consists of layers of an orthotropic material whose 
principal material axis makes an arbitrary angle e with the x axis (Fig. 
1). The constitutive equation (in condensed notation) can be written as 

= [C]'^({e*}'^-{a}'^AT-{6}‘^AH) (3.1) 

V I 

u 

[C] is the transformed 6x6 stiffness matrix 

ir 

(a) is a 6x1 stress vector 

l( 

{e*} is a 6x1 total strain vector 

Ir 

{o} is a 6x1 vector of thermal coefficients 
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Figure 1. fieornetry of Laminate with Idealized Damage Zone 
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{B> is a 6x1 vector of moisture coefficients 
aT is temperature change (from stress-free temperature) 
aN is moisture concentration change (from OX) 
and k refers to the kth layer. Explicit expressions for these quantities 
can be found in Appendix A. 

Noting that strains (and stresses) are functions of y and z 
only, it can be shown that suitable integration of the strain-displace- 
ment relations 

* **»x* V “ ^’y* V '**z 

Yy 2 ■ ^*2 ^*y* ^xz ~ * ^*x* ^xy ~ “*y * 

yields the following displacement field 

u - (C'jy ^ 

.2 

V = (C^z + Cj)x - 2 “ + V(y.z) (3.3) 

w = -C^xy + CgX -C 2 + H(y,2) 

The imposition of the proper symmetry conditions for axial loading 

u(x,y,z) = u(x,y,-z) 
v(x,y,z) = v(x,y,-z) 

w(x,y,z) = -w(x,y,-z) (3.4a) 

v(x,y,z) = -v(x,-y,z) 
w(x,y,z) * w(x,-y.z) 

as well as the experimentally determined condition (Pipes and Daniels 
[13]) 



8 


u(0,y,z) = -u(0,-y,z) (3.4b) 

results in the following displacement field: 

u = + U(y.z) 

V = V(y,z) (3.5) 

w = H(y,z) 


where = u,^^ = the uniform axial strain. In addition, the dis- 
placement boundary conditions on the midplane and centerline are 


U. 2 (y.O) = 0 U(O.z) = 0 

V. ^(y,0) = 0 V(0,z) = 0 (3.6) 

W(y.O) = 0 W.y(O.z) = 0 

This displacement field and the boundary conditions were first presented 
by Pipes and Pagano [14]. 

For completeness, the simplified equilibrium equations, with 
stresses independent of x and neglecting body forces, are shown below. 




ay 


+ 



= 0 


+ 

ay 


92 


(3.7) 


9t 9a 

9y 32 


0 


In addition, the full governing differential equations for the kth 
layer, as developed in reference [12], are presented below. 
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* Ss^’zz * W’yy * W*zz * ^^36 * ^45^***yz^*' “ ° 
^^26^*yy * ^45^ zz ^ ^22'^*yy ^ W*zz * ^^23 * ^44^“*yz^^ * ° 
^^^45 ^ Se^^’yz * ^^44 * ^23^^*yz * ^44**’yy * ^33‘^’zz^'' * ° 

3.2 Finite Elanent Formulation 

As in the previous section, only the main points of the formula- 
tion will be presented tore. A complete presentation can be found in 
references [1,2]. 

The finite element solution involves the subdivision (or discre- 
tization) of the structure into a finite number of elements (Fig. 2). 
For each elc ant, a set of interpo'ztion polynom'**ls is chosen to 
represent the displacement field within that element as functions 
of the element's nodal displacements. Using the strain-displacenent 
relations (3.2) and the constitutive relationship (3.1), the stresses 
and strains in each element can be found as functions of the nodal 
displacements, /^plying a variational principle, such as minimum 
potential energy, a set of equations relating nodal displacements to 
nodal forces can be obtained for each element in the form 

(3.9) 

where 

is the elemental stiffness matrix 

(a) 

(u)^ ' are the elemental nodal displacenents 
( 

{F}' are the elemental nodal forces 
These elemental relationships are then assembled into a global set 




o 


Y 


FSnure 2. Discretization of fiuarter Cross-Section of Damaged 
Fou*" layer Symmetric Laminate. 
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of equation*: relating all nodal displacements to all nodal forces. The 
application of boundary conditions and solution of these equations 
gives the displacements and thus the stresses and strains over the 
entire structure. 

Constant strain, constant stress triangular elements are utilized 
by NONCOM2. The interpolation polynomials have the form 

u * a^ + a 2 y + OjZ + 

V = a^ + a^y + a^z (3.10) 

w = a^ + agy + agZ 

When substituted Into the strain-displacement relations (3.2), these 
Interpolation polynomials yield constant strains over each elanent. 

The uniform normal strain Is ^ and the constants a^ through ag are 
functions of the element's nodal displacements and nodal spatial 
coordinates. 

Substitution of (3.10) into (3.2) yields the following strain- 
displacement relationship 

af^+cv2+evg 
bw^ +dW2+gWg 

bv^ +d V2+gv2+aw.| +cw2+ewg 
bu^+du2+gug 
au^+cu-j+eu^ 

where 

ORIGINAL PAGE IS 
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A = area of element 

UpU 2 >u^ = x-displacement at nodes 1,2,3 respectively 
Vi,V 2 ,V 3 = y-displacement at nodes 1,2,3 respectively 
WpW 2 ,W 2 * z-displacement at nodes 1,2,3 respectively 
and a,b,c,d,e,g are known constants Involving nodal spatial coordin- 
ates. 

For the case of uniform thermal load the strains are 



(3.14) 
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Ug = internal strain energy 


arJ 

Wg = potential energy of the applied loads. 

The internal strain energy for an element is 

Ue = ^J{€}^[C]{e)dVol (3.15) 

where [C] is as defined in (3.1) and {el is as defined 3.11) or 
(3.13), depending on whether it is mechanical or therma jroscopic) 
loading, respectively. For a unit thickness element with constant 
strains, reduces to 

Ug * ^ {e}^[C]t€) (3.16) 

and becomes 
e 

Wg = (3.17) 

where and are as defined in (3.9), 

Minimization of equ. (3.14) with respect to no ’al displacements 
yields the elemental stiffness matrix plus strain, thermal and 
hygroscopic related vectors. Exact forms for these can be found in 
reference [3]. The finite element meshes used for this report were 
generated by a program developed by Bergner, Davis and Herakovich 
[15] and can be found in Appendix B. 


3.3 Failure Analysis 


3.3.1 Tsai-Wu Failure Criterion 
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The three dimensional stresses: from N0NCCM2 provide a unique op- 
portunity to study a three dimensional failure criter-ion. Tsai and Wu 
[16] proposed such a criterion, postulating a failure surface of the 
form 


where the contracted notation is used and i,J=l,2,,..6. The strength 
tensors, F. and F... are of the second and fourth rank, respectively, 

» i J 

and F . . is assumed to be symmetric. 

* 0 

To determine F. and F^.^ in terms of uniaxial str*»ngths, a series 
of "thought" experiments are conducted. Assuming only a longitudinal 
stress on a unidirectional laminate, the tensile (X^) and compressive 
(X^) strengths are substituted into equ. (3.18), resulting in the 
following 


*"l^t '^ll’^t ^ ^ 


F,X + F,,X^ = 1 
1 c lie 


(3.19) 


Solving tt ‘‘e two simultaneous equations and F^-j are found to be 

(3 20) 


t C 


■1 


•11 

In d similar fashion, the following are found 

F = 1_. 1- 

^ ''t ''c 


- _ _ll_ 

22 Y Y 
c c 
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(3.21) 




1 




-1 


For the shear it is noted that shear strengths are inde- 

pendent of direction (and thus of sign) in orthotropic oaterials. Thus 
all terms involving s^iear to the first poNer are dropped. The only 
non-zero terms involving shear are then Fjg, and Fg^. With 
thought experiments similar to those used previously, the non-zero 
shear terms are found to be 



1 

(^^ 

1 


^^13^ 

1 

(s 


(3.22) 


where S-. refers to the shear strength in the i-j plane. 

The only terms left to be determined are F^^* ^23* 

ever, no simple uniaxial test will determine these interaction terms. 
They must be experimentally found by some bia.%ial test of strength. 

Wu [17] describes how to ''*^*ain optimum ratios of biaxial stress. Pipes 
and Cole [18] found, however, that vari?*'’ons 400 percent in the 
value of f ? interaction te m rr y »ittle difference in 

the 'trength envelope. 

For an orthotropic material, ther , the Tsai-Wu criterion takes the 
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form: 


F^O] + + F3O3 + F^^cJ ♦ F^^o^ 


+ F,,C3 + F^^T23 + f"55T^3 + ^66^12 


33 


(3.23) 


+ 2F^2°i°2 * ^^13°1°3 * ^^23°2°3 ^ ^ 


Since the Tsai-Wu criterion has a tensorial form, it is valid in all 
coordinate systems and the transformation has a simple tensorial form. 
For a 0 rotation about the 3 (z) axis, the Tsai-Wu criterion takes the 
form: 


* ’"2"/ * ''3'’2 * 

* ^22”l * fy-z * f«Vz * fs5'L 

* ^66'ly * ^""ie’x’xy * ^^26Vxy 


(3.24) 


* ^''M‘’z"xy * ^^ks'yi\z * ^^2\'’y 


* 2fl3'’x'’y * ° ’ 


The exact form of the F' terms can be found in Appendix C. 


3.3.2 Failure Model 

Using the Tsai-Wu criterion, the stresses on an element are used as 
input into equ's, (3.24). If the function has a value greater than or 
equal to one, the element has failed. However, a failed element is 
difficult to handle in a finite element solution. 

Sandhu [19] proposed that when failure occurs, the moduli of the 
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material should be set to large negative numbers. This essentially 
"unloads" the failed material . The major drawback with this procedure 
is the creation of negative diagonals in the global stiffness matrix, 
causing an unstable matrix. The method used in references [1-3] is to 
reduce the moduli by the percent of overstrain in the failed direction, 
a ccmplicated procedure and one that still allows a failed element to 
carry load. 

A simpler approach, and the one adopted in this study, is to 
reduce all the moduli in a failed element to some small value. In 
addition, the stresses this element are set equal to zero. By 
reducing the moduli and setting the stresses equal to zero, the 
failed element carries no load and does not contribute to the overall 
stiffness of the laminate during subsequent loading. However, even 
this approach can cause an unstable stiffness matrix, as will be 
shown later (see Section 5.4.2). 



Chapter 4 

THE FINITE ELEMENT PROGRAM (NONCOM2) 

The program NONCOMl , a modified version of NONCOM, was modified 
further for this study and called N0NC0M2. Changes made for this 
version include a more efficient equation solver, a failure criterion 
based on the Tsai-Hu function (Chap. 3) and an improved handling of 
hygro thermal effects. 

4.1 Equation Solver 

The program NWICOMI used an equation solver called SEESOL [20], 
which stored the stiffness matrix in blocks on tape. However, I/O 
charges on VPIiSU's IBM 370 make this an expensive method of solving 
simultaneous equations. Since the VPI4SU system has an extremely large 
core and can thus store large matrices, a more economical inethod is an 
in-core solution. 

An equation solver, COLSOL, presented by Bathe and Wilson [21], is 
an entirely in-core soluLicn. The upper triangular half of the sym- 
metric stiffness matrix is stored, colunm by colunm, in a row vector, 

A. For a banded matrix, terms above the "top" of the bandwidth are 
zero. In order to minimize storage requirements, terms above the "top" 
of a column are not stored. This is an extremely efficient method of 
storage for symmetric, highly banded matrices similar to the type shown 
below. 
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After the global stiffness matrix is stored, Cholesky decomposition 
is used to solve the equations and provicj the nodal displacements. 

4.2 Material Considerations 

The program N(WC0M1 had several restrictions on the temperature and 
moisture dependence of material properties. Only elastic and shear 
moduli, and thermal and moisture coefficients were assumed temperature 
and moisture dependent; all other material properties were assumed con- 
stant. Also, it was assumed that tensile and compressive moduli varied 
in the same manner. The modified version, NONCOM2, used in this study 
allows all material properties to vary as functions of temperature and 
moisture. Further, tension and compression values are allowed to vary 
as independent functions. 
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Hygrothermal data is presently found in the literature in the form 
shown in Fig. 3. In the program, each of the curves is modeled as a 
series of linear segmented lines. Values of percent retention of a 
material property are found by linear interpolation between the end- 
points of the segments. 

The material system used in this study was T300/5208 graphite-epoxy. 
Complete mechanical and thermal data were found in a report by Hofer, 
Larsen and Humphreys [22], Hygroscopic data for this material can be 
found in reference [22] and a report by Kiebler [23], but the data are 
neither complete nor consistent. Since the loadings for this study were 
all mechanical only, the results were not affected by the incomplete 
hygroscopic data. 

Material property values can be found in Appendix D. In addition, 
a user's guide for the input of hygrothermal properties into N0NCCM2 can 
be found in /^pendix E. 

4.3 Nonlinear Analysis 

The finite element analysis of Chapter 3 requires linear elastic 
properties. Since composites often have nonlinear mechanical proper- 
ties, an incremental procedure is utilized to account for material 
nonlinearities within the framework of a linear elastic analysis. 

4.3.1 Incremental Procedure 

Using the incremental procedure, the mechanical properties are 
linearized for each increment and the load, either thermal, hygroscopic 
or uniform axial strain, is applied as a series of increments. Sum- 




TEMPERATURE 


Figure 3. Typical Percent Retention Curves as a 
Function of Temperature and Moisture 
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(nation of the linear solutions of stresses, strains and d1sp'ocen®r.ts 
for each increment yields the total solution. At each step of the 
procedure, the mechanical properties of each element are altered de- 
pending on the temperature, moisture and/or strain level. This prx)- 
cedure allows the properties of an element to vary independently of all 
other elements. 

4.3.2 Determination of Tangent Modulus 

Nonlinear stress-strain curves were represented by Ramberg-Osgood 
[24] approximations of the form 

n. 

e = ^ + K^a ^ i = 1 or 2 (4.2) 

where E is the elastic modulus and K- and n. ar^ Ramberg-Osgood coef- 
ficients. The procedure to calculate the four constants K. and n. can 
be found in reference [1]. The tangent modulus can be defined as 

K.En.o +1 

p 

Noting Fig. 4 the corresponding stress, a , at the end of load increment 
P is 

P 

= r (4.4) 

where Ae^ is the increment of strain and the tangent modulus during 
the jth load increment. For the P+lth increment, the tangent modulus 
becomes 
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Fiqure 4. Determination of Tangent Modulus with 
Rambero'Oscood Approximations 
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= p i = 1 or 2 (4.5) 

K,En,[ V ArjEj]"i’’'+1 
' ' j=l 

Using equation (4.5) and principal material strains, the tangent modulus 
for the next increment is calculated at the end of each increment. Fnm 

p 

Fig. 4, it can be seen that the strain at which the tangent 

modulus, is calculated, differs from the strain actually 

applied. This difference can be made arbitrarily small by choosing a 
sufficiently small increment. For the computer analysis, equations of 
the form (4.5) were used to determine the tangent moduli ^ 22 ’ ^ 33 ’ 
^23’ ^13* ^12 finite element. It is assumed that shear 

response is independent of sign while extensional behavior can be 
different in tension and compression. 

4.4 General Notes 

In Appendix D, the values of the F^ 2 » G ]3 *^23 inter- 

•10 2 

action terms are given as -0.58x10 /(PSI) , a value taken from 

-10 2 

reference [18]. Wu [17] suggests a value of 2.0x10 /(PSI) as more 
appropriate for graphite-epoxy. Several computer runs showed, however, 
that the differences in the two values did not affect the overall 
trends. Since reference [17] was found after all the computer runs 
had been made the cases presented in Chapter 5 were not rerun with Wu's 
value. 

The initial tangent modulus for all three versions of NONCOM is 
chosen as the tensile elastic modulus where, for compressive loading, 
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compressive moduli would be more appropriate. However, this factor did 
not significantly alter the overall results or conclusions. 



Chapter 5 

RESULTS AND DISCUSSION 

The analysis presented in Chapters 3 and 4 was used to study 
various laminates of differing ply orientations. Four and eight layer 
laminates with or without cutouts were analyzed. Because of symmetry 
only a quarter of the laminate cross-section, shown in Fig. lb, was 
examined. Both linear elastic and non-linear results are p»*esented. 

5.1 Averaging of Finite Element Results 

The finite element analysis of Chapter 3 is b?sed on a displacement 
formulation with linear interpolation functions over each element. This 
lormulation results in stresses and strains which are constant over each 
element. Constant stresses for an individual element give a distri- 
bution of stresses which appear discontinuous over a series of elements. 
In order to produce smooth distributions, stresses are averaged over 
several elements. 

Noting Fig. 5a, stresses for a through- the- thickness distribution 
are averaged along the line A-A. The stresses presented at point E 
correspond to an average of the stresses in elements 1 and 2. This 
method of averagim was used in all through- thickness plots. 

For interfacial distributions, the a , t and t stresses ar*e 

c jL XZ 

averaged along the interface between the two layers, shown as line C-C 
in Fig. 5b. Since these three stress components must be continuous 
across the interface, the results are an average of elements above and 
below the interface. Thus, for pcint H the stresses are an average of 
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(a) THROUGH -THICKNESS 
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Figure b. Averaginq of Finite Element Results 
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elements 14, lb, 16 and ^7. The stress components a , c and t .which 

X y x> 

are not necessarily continuous across an interface, are averaged along 
line B-B and line D-D. At point F, stresses are averaged over elements 
8 and l while at point 6 stresses are averaged over elements 10 and 11. 

5.2 Convergence of Finite Element Results 

The ^inite element analysis, as mentioned previously, is based on 
a displacement formulation. However, the boundary conditions at the 
cutout and free edge are stress-free boundary conditions. Due to the 
approximation introduced by the displacement formulation, as well as 
the use of constant stress, constant strain elements, the stress-free 
boundary conditions will no*^, in general, be satisfied. However, 
Herakovich [25] has shown that the free edge conditions „ond to be met 
as smaller and smaller elements are used near i^Se free edge. Thus, 
convergence of the finite element solution al a particular pcint is a 
function of mesh size and displacement gradient. 

Satisfaction of the equilibrium conditions for sections of the 
laminate (Fig. 6) is also a function of mesh size and displacement 
gradients. Similar to the stress boundary conditions at a point, 
equilibrium conditions for sections will tend to be met as the size of 
the elements is reduced. Equilibrium for the body as a whole, though, 
is automatically satisfied due to the formulation. 

Thus, no finite element mesh based on a displacenient formulation 
will satisfy stress boundary conditions or equilibrium conditions 
exactly. However, by using meshes with similar sized elements at cutout 
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Figure 6. Partial Free Body Diagram of Quarter Cross-Section 
of Damaged Four Layer Symmetric Laminate 
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and free edges (areas of high displacement gradients) a sufficiently 
adequate stress distribution fo»* relative comparisons and general con- 
clusions is obtained. 

5.3 Linear Elastic Results 

This section contains linear elastic results for laminates with 
aspect ratios b/H = 68.2 (four layer laminate) and b/H = 34.1 (eight 
layer laminate) (Fig. 1) and applied uniform axial strain ^ = -0.1 
percent. ITie curves in this section were drawn by the VPI4SU CALCOMP 
plotter. 

5.3.1 Unidirectional Laminates 

The stress state in unidirectional laminates with all fibers 
oriented at an angle e is a uniform uniaxial stress, o^, with or without 
a cutout. The value of is equal to the laminate value while the 
other stress components are approximately zero. Thus, for unidirec- 
tional laminates, the cutout does not alter the uniform stress state and 
produces no stress concentrations. 

5.3.2 Four Layer Laminates 

5.3.2. 1 The [0/90]^ Laminate 

The Oy stresses for the [0/90]^ laminate with one layer removed as 
the idealized damage zone are shown in Fig. 7. These stresses ara shown 
along the interface from the cutout edge (point A) to the free edge 
(point B). The stresses for the same layup with no layer removed are 
shown by the dashed lines. Near the free edge, the variation of as 
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Fiqure 7, along Interface of [0/90]^ 
Laminate {^=-0.1") 


(MPa) 
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well as all other stress components, is exactly the same with or without 
a cutout. Also, because of the dimensions of the specimen, the 
o (and o and t ) stresses in the interior region (away from the 

Jr A AJr 

edges, between the cutout and free edges) are constant and equal to the 
theoretical laminate values for [0/90]^ laminate (Table 1). 

The and stresses, along the interface, shown in Fig. 8 and 
9, respectively must satisfy several equlibrium conditions. T.ie partial 
free-body diagram (Fig. 6) shows the 0® (Cj=0“) layer of a [0/90] 
laminate with a section removed and a cutting plane, C-C, passing 
through the interior region. In order to satisfy force equilibrium, the 
follow'*! be true. 

f ' 0 



T _dy 
yz 


0 



o^dy = 0 




o^dy = 0 



(5.1) 
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lABLE 1 

Comparison of Lamination Theory and Finite Element 
Results in Interior Regions (y/b=0.67) for [0]/®2^s 



(psi) 

(psi) 


®1 

®2 

®1 

«2 

®1 

«2 

[0/90]j 

Lamination 

-27790 

-1552 

-535 

535 

0 

B 

N0NCGH2 

-27760 

-1551 

-535 

535 

0 

0 

[±4S]j 

Lamination 

-3703 

-3703 

0 

B 

-1588 

1588 

NONCOM2 

-3703 

-3703 

0 

0 

-1588 

1588 
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Oy is the theoretical laminate value, 
the following must also be true. 


To satisfy moment equilibrium. 


jf oprdy = a‘-y(h^/2) 
y® ci^dy = -Oy{h^/2) 


(5.2) 


For the reasons cited in Section 5.2 and because the stresses are not 
known at the exact edge, these conditions are not well satisfied. 
However, by assuming a linear stress distribution from the last point 
given by the finite element method to the edge of the laminate, the 
equilibrium conditions 5.1 and 5.2 can be satisfied within two percent. 
A larger number of elements at the edges would also tend to satisfy the 
conditions more accurately. 

Each of the Figs. 7, 8 and 9 show the free and cutout edge effects 
extending approximately 14 h^ into the interior of the laminate. The 
free edge stresses are slightly larger than the stress values at the 
cutout edge. Thus, for this laminate, the free edge is more critical 
than the cutout edge. 

The value of the Tsai-Wu function is largest at the free edge, 
again showing that the free edge is most critical. The through-the- 
thickness variation of the Tsai-Wu function near the free edge. Fig. 10, 
indicates that the largest value occurs near the top (z/hQ= 2.0) of the 
laminate. (The values of the Tsai-Wu function in Fig. 10 have been 
scaled to be between 0 and 1.) Examining the terms of the Tsai-Wu 
function for a 0” layer (Equ. 3.23), the critical term was determined to 
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Through-Thickness Distribution of Tsai-Wu Function 
for [0/90]^ Laminate ( y/b=0.99) 


Figure 10. 



be f 2 ^y(Oy= 02 )’ Noting Fig. 11. which shows the through-the-thickness 
variation of the critical term will have the largest absolute value 
near the interface. Since in the 0” layer is negative, F 2 a^ is less 
at the interface than near the top. This indicates failure may begin 
near the free edge and near the top of the [0/90]^ laminate. 

5. 3. 2. 2 The [*45]^ Laminate 

A cutout in the [-45]^ laminate produces a number of important 

effects, particularly near the cutout. Figure 12 shows in the +45" 

layer rising to a value 124 percent above the stress in the interior 

region. Fran a value of zero in the interior, the stress (Fig. 13) 

attains a value of approximately -1000 psi near the cutout. The dashed 

line in Fig. 13 indicates the variation of in the -45® layer fron 

point E 1.0 point The r stress (Fig. 14) in the +45® layer near the 

xy 

cutout rises to a value 160 percent above the lamination value in the 

interior. The interior region values for the c , c and r stresses, 

X y xy 

which are the same as the theoretical lamination values, are given in 
Table 1. Again, as in the [0/90]^ laminate, free edge distributions for 
all stresses are the same with or without a cutout. 

The c, stress. Fig. 15, shows an extremely important effect of the 

4 . 

cutout. Moment equilibrium (Equ's. 5.2) requires that 
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since o'" is zero in the interior reqion. The o variation near the 
y ^ z 

cutout shOMS the correct trend, having three envelopes to produce a zero 


moment. The variation near the free edge, however, does not show the 
correct trend due to a steeper displacement gradient in this area. In 
addition, the variation shows the effect of the cutout edge extending 
twice as far into the interior region as the free edge effect. 


The and shear stresses, shown in Figs. 16 and 17, respec- 
tively, also exhibit the influence of the cutout edge. The stress 
near the cutout is very large in comparison to the same stress near the 
free edge while the stress near the cutout is only slightly larger 
than the free edge value. In addition, the distribution has two 
envelopes near the cutout whereas the single envelope at the free edge 
is the distribution expected. Hot-.tver, force equilibrium of the left 


top section of Fig. 6 between cutting planes JJ and CC requires that 


/; -r -i: 


x„dY - 0 


(5.4) 


Since t is greater along line JJ than CC (Fig. 14) and by convention 
xy 

in the opposite direction, the integral of must be negative (and 
thus, by convention, in the positive direction). As the cutting plane 
JJ is ntoved closer to he cutout edge, \ goes to zero and the integral 

of ; becomes positive to balance the t on plane CC (a similar • 

xy Xi 

stress distribution should exist near the free edge but doesn't due to 
the coarseness of the grid). Similar arguments can be advanced to 
explain the two envelopes of the distribution (Fig. 16) involving 
the o stresses. 

y 
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The Tsai-Hu function along the interface (Fig. 18) indicates the 
largest value occurs in the interior of the -45® layer near y/b ■ 0.4. 
The through-the-thickness variation near this point (Fig. 19) shows the 
largest value occurring not at the «45 interface but near the midplane 
(z/h^ * 0.0). The through-the-thickness variation of stresses Indicates 
that the stress (Fig. 20) is the principal cause of the maxiiaum Tsai- 
Wu va1»e occurring at the midplane. The stress rises to a maxinun at 

the midplane in order to balance the moment produced by the o stress 
(Fig. 21) along the plane y/b - 0.4. Thus, failure may initiate at the 
midplane and a small distance in from the cutout. 

5.3.3 Eight Layer Laminates 

Partial free body diagrams of the quarter section of eight l^r 
laminates with cutouts are shown in Figs. 22 and 23. There are three 
interfaces, referred to in the following manner; the upper interface 
from point A to B, the middle interface from point C to 0, and the lower 
interface from point E to F. Also, the interior region refers to the 
section away from the free and cutout edges where the line KK is shown. 
In the interior region, the laminate solution for a 
laminate holds. The laminate solution also holds ne^r the centerline 
(y=0), but for a depending on whether 

one or two layers is removed. 

5.3.3. 1 The [02/^45]^ Laminate 

The Oy stress distributions for a one layer cutout and for a two 
layer cutout are presented in Figs. 24 and 25 respectively. Each shows 
the value of changing from the centerline laminate value to the 
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Throuoh-Thickness Distribution of Tsai-Wu 
Function for [*45}^ Laminate ( • y/b=0.40) 


Figure 19. 
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Figure 20. Throuah-Thickness Distribution of for a [-45]^ 
Laminate (f^=-0.1", y/b=0.40) 
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interior region laminate value. The finite element results in those 

regions for o and x , as well as a , are exactly equal to the laminate 
X xy y 

values. The presence of the cutout, whether one or two layers, does not 
significantly increase the mag»'itude of those stress components at the 
cutout edge. Further, all stress components at the free edge have the 
same variation \ th or w.thout a cutout. 

The cutout does have an unusual effect on the distributions as 
shown in Fig. 26, which gives the results for a laminate with a two 
layer cutout. Along the lower interface a couple is produced to balance 
the moment produced by the Oy stresses in the interior region (line 
KK, Fig. 23). Though not as high as the stresses piouuced at the free 
edge, this resultant o, couple in the material just below the cutout 
edge (y/b = 0.33) could initiate Tailure if a flaw existed in this 
region. 

Typical x^^ and x^^^ distributions are shown in Figs. 27 and 28, 

respectively, for a laminate with a one layer cutout. The distribution 

for a two layer cutout are similar except that (as expected) th® middle 

interface distribution resembles the upper interface distribution. For 

the Xy^ distributions (Fig. 27), the envelopes near y/b = 0.33 on the 

middle and lower interfaces are necessary for force equilibrium (eFy = 

0). For x^^ ( ig. 28), the double envelopes near the free edge on both 

the upper and middle interface are caused by a slight rise in t near 

xy 

the free edge (by reasoning similar to Section 5. 3. 2. 2). 

5. 3, 3. 2 The [+45/02]^ Laminate 

The stresses, shown in Figs. 29 and 30 for one and two layer 
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Figure 26. along the Interfaces of Laminate with 

Two Layer Cutout ( f^=-0.1'^). 
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cutouts are the most Important stresses for this laminate. As can be 
seen from both figures, the tensile stresses along the interfaces 
below the cutout are the highest tensile stress encountered. In fact, 
the tensile stress on the lower interface of the [ 145 / 02 ]^ laminate 
with a two layer cutout (Fig. 30) is the largest tensile stress 
near the cutout for either a [*45/0,]. or [0./t45]_ laminate with or 
without a cutout. As before, this suggests failure may occur along this 
interface, particularly if it is a flawed or damaged region. 

The stress distributions along the interfaces for one and two 
layer cutouts (Figs. 31 and 32, respectively) have several interesting 
differences. Along the lower interface, the envelope just below the 
cutout peaks at a higher value in the one layor cutout than in the two 
layer cutout. Further, the negative envelope near the cutout along the 
upper interface of the one layer cutout is not present in the distri- 
bution along the same interface in the two layer cutout. A similar 
pattern occurs in the distribution (Fig. 33) along the upper inte- 
rface where the tvio envelopes near one layer cutout do not occur in the 
two layer cutout. Again, arguments similar to those presented in Section 

5.3. 2.2 can explains these patterns, particularly since t and 

xy 

Oy near the cutout in the ♦45” layer rise to values 142 percent and .?75 
percent, respectively, above the interior region values for a laminate 
with a one layer cutout. 

The Tsai-Wu function for this laminate attains its maximum value in 
the zero degree layer near y/b * 0.4 for both one and two layer cutouts 
as shown in Figs. 34 and 35. This point concides with the point of 
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maximum stress (Tigs. 30 and 31) and, when the terns of the Tsai-Wu 
function are examined, it is found that i.he stresses are the major 
contributing factor. Thus, failure may initiate at an internal point 
below the cutout and away from tne free and cutout euges. 

5.3.4 Quasi -Isotropic Laminates 

The quasi -isotropic laminates studied were [0/*45/90]^, [90/i45/0]^ 
and [i45/0/90]^, all eight layer laminates. The commeiics pertaining 
to interfaces and interior regions made in the introductic to Section 
5.3,3 are, of course, valid for these special cases of eight layer cutouts 
(Figs. 22 and 23) were examined for all three laminates. 

5.3.4. 1 The [0/t45/90]^ Laminate 

The a (Fig. 36), a (Fig. 37) and x (Fig. 38) stress distribu- 
X y xy 

tions for a laminate with a two layer cutout indicate that the cutout 
produces more severe stress concentrations when it is along the ±45 
interface. In the +45 layer, and near the cutout, stress concentra- 
tions of 1.2 for 0 , 1.34 for o and 1.31 for x are present. As be- 
X y xy 

fore, the interior region values are the lamination theory values for 

a [0/+45/90]^ laminate. It should also be noted that the trends for 

0 and X at the free and cutout edges are consistent with the 
y xy 

boundary conditions, as both stresses tend to zero near the edges. The 
dashed lines in the figures indicate the differences for a one layer 
cutout. 

The o^, and x^^ interlaminar stress distributions for the two 
layer cutout are shown in Figs. 39, 40 and 41, respectively. The 
maximum positive occurs along the lower interface near the free 
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e<l9^. but the value below the cutout edge Is nearly equal to this max1> 
mum value (Fig. 39). Equilibrium considerations can explain the double 
envelopes of the and distributions (Figs. 40 and 41) near the 
cutout as In the (^45]^ laminate (Section 5. 3. 2-2). For one layer cutout 
gradient distributions are not present near the cutout due to the 
coarseness of the mesh (see Section 5.2). Free edge stress distributions 
are the same, though, for one and two layer cutouts. 

The distributions of the Tsal-Uu function along the Interface of 
a laminate with a two ' yer cutout (Fig. 42) Indicate that failure Is 
predicted to occur In the 0* layer just above the 0/45 Interface near 
the free edge. (The distributions for a one layer cutout are shown by 
the dashed lines.) Examination of the Individual terms of the Tsai- 
Vki function indicates that no one stress component Is the critical 
term. With the exception of all stresses exhibit free edge values 
which contribute to the maxiimmi Tsai -Mi value occurring in the 0* 
layer near the free edge. 

5. 3. 4. 2 The [90/*45/0]^ Laminate 

The major difference in stress distributions between this stacking 
sequence and the [0/*45/90]^ laminate is shown by the distributions 
(Fig. 43) for a two layer cutout. Fig. 43, which is for the cutout 
a’ong the *45 interface, shows in the *45* layer tending to zero near 
the cutout. In contrast. Fig. 37 exhibits an increase in the magnitude 
of in the M5" layer near the cutout for the [0/*45/90]^ laminate. 

The interlaminar shear stress is also affected by the stacking 
sequence, as evidenced by a comparison of the middle interface stress 
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Figure 43. along the Interfaces of [90/*45/0]l^ Laminate 
with Two Layer Cutout (f «-0.1*l. 



(MPo) 





distributions of Figs. 40 and 44. 

The Oy stress distribution for a one layer cutout (Fig. 45) also 
shows a difference between this stacking sequence and the previous 
one. As indicated in Fig. 45, in the top (90'') layer increases 
120 percent above the interior region value. This is i’l contrast to 
Fig. 37, which does not show this large increase for the [0/ *45/90]^ 
laminate. However, the corresponding interlaminar stress (not 
shown) for a [90/ *45/0]^ laminate is not affected significantly by 
this increase. The i , distribution for the upper interface does not 
exhibit a double envelope near the cutout edge, but rather a distribu- 
tion similar to that near the free edge. 

Although no* shown in the figures, the t stress distributions for 
a [90/ *45/0] laminate with a one or two layer cutout follow the same 
pattern as the distribution shown in Fiq. 44. The interlaminar 
shear component (Fig. 46) is affected in the same manner as was the 
component. The double envelope near the cutout, seen in Fig. 41 
along the middle interface, is not present in Fiq. 46. Thus, the cuiout 
does not act as a stress concentration for ; and produces no unusual 

Ajr 

effects in **■ distribution. 

Fig. 47 presents results for the interlaminar nctinal stress 
for a laminate with a two layer cutout. In contrast to the results 
for a [0/'4S/9il]^ laminate (Fig. 391. the stresse^ along the upper 
interface are larqe for the laminate. Also, the 

stresses near the edges are positive for the [90/t4S,'0]^ laminate and 
negative for the [0/'4o/go]^ laminate. Thus, the [*50 '*45/0]^ laminate 
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Is more susceptible to the delamiudtion mode of failure at the edges. 

The distributions of the Tsai-Wu function for the [90/±45/0]^ 
laminate with one and two layer cutouts (Figs. 48 and 49, respect^'vely) 
indicate that failure is predicted to initiate near the free edge in the 
0® layer. However , it should be noted that high values of the Tsai-Wu 
function are also present below the cutout in the 0® layer of the 
laminate with a one layer cutout (Fig. 48) and near the cutout in the 
+45® layer for the two layer cutout (Fig. '•*'). Thus failure could 
initiate at the free edge in the 0® layer and, at a higher strain level, 
in either of the two regions near the cutout where high values of the 
Tsai-Wu function occur. The nonlinear results to be discussed in 
Section 5.4.2 show this behavior. 

5. 3. 4. 3 The [±45/0/90]^ Laminate 

The [±45/0/90]^ laminate provides little new insight. Like the 

other two quasi-isotropic laminates, the distributions for a two 

layer cutout (Fig. 50) show approaching the zero boundary value 

rapidly near an edge. Although not shown, a similar trend occurs in the 

T distribution. Honkcver, when the bottom of the cutout is along the 
xy 

±45 interface \ me layer cutout) o and t exhibit stress concen- 

y 

tration factors of 1.27 and 1.26, respectively, near the cutout and in 
the +45° layer. As before, these stresses induce a double envelope near 
the cutout for the x and x distributions along the upper interface 

y2 X . 

(also not shown). 

5.4 Inelastic Results 

The [90/±45/0]^ laminate, with and without cutouts, was stuoied for 
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inelastic behavior. The axial compressive strain was applied incre- 
mentally to failure in increments a ^ = -0.05 percent strain. Failure 
for the laminate was defined as the failure of all elements using the 
Tsai-Wu criterion for individual elements (Section 3.3.2). In the case 
of the one layer cutout, the stiffness matrix became unstable and the 
solution procedure failed. However, since the load had already dropped 
appreciably, chis was taken as an indication of total failure. 

5.4.1 Average Stress-Strain Results 

The stress-strain behavior up to the first large decrease in stress 
is shown for all three laminate configurations in Fig. 51, where is 
the average compressive stress over the entire cross-section. As Fig. 

51 indicates, the laminate modulus is approximately the same for the 
one layer cutout and the complete laminate, but is higher for the two 
layer cutout. This is due to the higher volume fraction of 0“ material 
in the laminate with the two layer cutout. Prior to failure, the 
modulus of each curve falls off slightly. In each case, this reduction 
in modulus corresponds to the failure of one or more large elements in 
the 0° layer in the interior of the laminate away from the edges. 

The total stress-strain curves are shown in Fig. 52. Though the 
maximum stress differs the maximum load (in units of force) for both 
one and two layer cutouts is 60 percent of the maximum load for the 
complete laminate. Also, for ^ between 0 and -1.4 percent, the various 
reductions in correspond mainly to failure of elements in the 0° 
layer. In all cases the 0“ layer failed before the other layers failed. 

In the case of the one layer cutout, the elements with nodes along 
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the centerline failed at = -1.^ percent. Since boundary conditions 
are applied along the centerline and midplane, this results in an un- 
stable stiffness matrix. In the cases of the two layer cutout and the 
complete laminate, the elements with nodes along the centerline failed 
at ^ between -1.4 and -4.0 percent when only a few isolated elements 
remained. Thus, no unstable matrix resulted. 

5.4.2 Failure Progression 

In all three cases, the free edge in the 0® layer was the first 
region to fail. As the strain level was increased, elements in the ±45° 
layers near the free edge failed. At the strain level ^ = -0.4 per- 
cent, several eleinents in the 0° layer near the centerline failed and, 
as the strain level was increased further, elenrents between the free 
edge and centerline in the 0° layer failed. Finally, with the 0° layer 
failed and elements at the edges in other layers failed, the laminates 
were essentially failed. 

In the region near the cutout there was a significant difference in 
the order in which elements failed, depending on whether the cutout was 
one or two layers. For the two layer cutout, elements failed along the 
cutout edge in the +45° layer and then the 90° layer. The pattern was 
similar for the one layer cutout but progressive failure also occurred 
in the +45° layer around the cutout. First, elements near the cutout 
failed and then, as the strain level increased, more and more elements 
around the cutout failed. 


5.4.3 Failure Prediction 
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The first failure of the laminate with a two layer cutout occurred 
at ^ = -0.4 percent. The distributions of the Tsai-Wu function for ^ 

= -0.35 percent for both one and two layer cutouts are shown in Figs. 53 
and 54. These distributions are similar to the elastic results pre- 
sented in Figs. 48 and 49, where ^ = -0.1 percent. Thus, the elastic 
results for a single increment give excellent predictions of where 
failure will occur. Both Figs. 53 and 54 indicate fialure will initiate 
in the 0° layer near the free edge for all three laminate configura- 
tions. This is as expected since the free edge stress distributions are 
the same with or without cutouts. Initial failure actually occurred in 
this region on the next increment ( ^ = -0.4 percent). 

Results for a one layer cutout (Fig. 53) show high values of the 
Tsai-Wu function below the cutout in the 0°, +45® and -45® layers while 
the results for the two layer cutout (Fig. 54) show high Tsai-Wu values 
near the cutout in the +45 layer only. In both cases, subsequent 
failure occurs in these regions at higher strain levels. 












TSAI'WU FUNCTION 



00 02 04 06 08 lO 

y 

b 

riqurr ^4. Pi ■strihutioti of T<;ai-Wii Function alonq the Interfaces of 
[00/ -4^/0]^ Laminate with Two Layer Cutout ( t^=-0.35' ). 







Chapter 6 
CONCLUSIONS 

This study was concerned with linear and nonlinear stress and 
failure analysis of an idealized damage zone in a laminated composite 
under compressive loading. The following conclusions can be made from 
results of the analysis. 

1. A boundary layer exists around the damage zone; it is 
similar to the boundary layer along the free edges. 

2. Interlaminar stresses in the boundary layer around the 
damage zone are similar to those produced at the free 
edge. 

3. The interlaminar stresses around the damage zone may be 
as large or larger than those at the free edge depending 
upon the fiber orientations a d stacking sequence. 

4. The boundary layer associated with the damage zone extends 
as far or farther into the inte'*ior of the laminate as does 
the free edge effect, depending on fiber orientation 

and stacking sequence. 

5. Only laminated materials with more than one fiber 
orientation are affected by the damage zone. 

6. Lamination theory is valid in interior regions away 
from fre? edges and the damage zone. 

7. The presence of a damaged zone induces interlaminar 
stresses in an interior region away from free edges and 
the damage zone. 
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8. Due to interlaminar stresses produced in an interior 
region, failure, defined by a Tsai-Wu failure criteriv.i, 
may initiate in the interior away from free edges and 
the damage zone. 

9. The cutout extending down to a ±45 interface may pro- 
duce large stress concentrations near the edge of the 
damage zone, depending upon the fiber orientations 
and stacking sequence. 

10. Elastic results adequately predict where failure 
initiates for the laminate considered. 

11. A damaged [90/ ±45/0]^ laminate fails at approximately 
60 percent of the load carried by an undamaged laminate. 

12. For a [90/±45/0]^ laminate, initial failure occurs at 
the free edge in the 0° layer. Subsequent failure in 
the interior of the 0° layer as well as at the free 
and cutout edges of the other layers leads to total 
failure. 

This investigation indicates that future areas of study should 
include: 

1. Linear and nonlinear analysis of hygrothermal effects 
in laminates with idealized damage zones. 

2. A three layer damage zone (cutout) in an eight layer 
laminate. 

3. The effects of varying geometric parameters such as 
cutout wic '■» and layer thickness. 
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4. More detailed nonlinear analysis. 

5. An improved method for modeling failure and unloading 
of an element. 

6. Experimental study of specimens with ide.»l'. ed damage 
zones to compare with theoretical analysis. 

7. Analysis capability for bending of unsymmetric lam'nates. 
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APPENDIX A 

CONSTITUTIVE RELATIONS 

The constitutive relationship for an orthotropic material in the 
principal material directions is 

{oil = [C]({e}^-{o}^AT-{B}^Af1) 

where 
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For a e rotation about tne 3 (z) axis (Fig. 1), the constitutive 
relationship becomes 


where 


{o} = [C]({€}-{a}AT-{e>AM) 


c c 
hi h? 


0 

0 

^16 

c 

h2 

ha 

0 

0 

^26 


'33 

0 

0 

^36 



'44 

'45 

0 

Synwetrlc 



^55 

0 
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and the various matrix and vector terms as functions of the principal 
material properties are given below (m=cose, n=sine). 


= m^C^^+2m^n^(C^2^2Cgg)+n^22 
C^2 = m^n^(C^^+C22-4Cgg)+(m^+n^)C^2 
^13 ^ ^13^" ^23 

C]g = -mn[m^C^^-n^C22-(m^-n^)(C^2^^^66^^ 
^22 ^ n^C^^+2m^n^(C^2^2^56^^^^22 

C^j = n^C,3«.^C23 

^33 " ^33 

C36 = mn(C23-C^2) 

2 2 

^44 ^ *" ^4'^" ^55 
C45 = mn(C^^-Cgg) 

C55 = n^C,X<:55 

^66 ° '"^"^(Cll*C22-2C32)*(m^-n^) Cjg 
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= m Oi+n 02 
Oy = n o-j+m ©2 

“2 

“xy * *211111(0^-02) 

®x “ n»^6^+n^62 
2 2 

By = n^e^+m‘-B2 

8, =63 

Bxy = -2mn(B^-B2) 
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APPENDIX B 

FINITE ELEMENT MESHES 

The following are the various linite element meshes used in this 
study. Each models the quarter plane of symmetry (Fig. lb). Thus, a 
two layer mesh mooels a four layer symmetric laminate and a four layer 
mesh models an eight layer symmetric laminate. Simllarlly, a one layer 
cutout means one layer removed on both top and bottom of the laminate (2 
layers total removed). 
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\ Fi.i, Two l.^yer Finitf- riewrit fV?',*-! for Ou/^rtor Cro^', -‘’/Oot i on oF Four 

Layer ‘iyrnf-trir Laninate. 



Figure B.?. Two Layer Finite Element Mesh for Duarter Cross-Section of Four 
Layer Symmetric Laminate with One Layer Cutout. 
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Foor l.aye»" Finitf* Flr-ffir-nt Mfj'ih for Ouarttir Cro'. -Sf'C tion of 
fiqht lH/<>r 'yrmotric l.aninato. 
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Figure B.4. Four Layer Finite Elenent Mesh for Quarter Cross-Section of 
Eight Layer Symnetric Laminate with One Layer Cutout. 


z 



Fiqut'f' B.5. Four l.riyer Finit*'* Flenent Me'^h for Quarter Cross -Soot ion of 
Fioht layer Syrrinetric Laninato with Two Layer Cutout. 
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APPENDIX C 

TSAI-WU FAILURE CRITERIA 


■’’he Tsai-Wu failure criteria has the tensorial fom (for an 
orthotropic material in the principal material directions) of 


■^''3j®3'^‘’44^L^''55^13^’'66*12 


+2F 


12 ^ 1^2 


,+2F, 


13°l‘^3'^^’^23°2®3 ^ 


(C.l) 


where the F. and F^^ terms are as previously defined in Chapter 3. 
For a 6 transformation about the z (3) axis (Fig. 1), the Tsai-Wu 
criteria takes the form 


+F‘ t^+ 2F' o T +2F’ o T 
66 xy 16 X xy 26 y xy 

+?Fl^cr T **'2F1 cT t +2F^«a a 
j 6 z xy 45 yz X2 12 x y 


* 2 Fj 30 ^ 0 y* 2 F' 3 OyO^ ‘ ' 


(C.2) 


whc.'e the F' terms, as f^iictions cf the unprimed F's and e, are as 
follows ("'. = cose, n = sir^*) 

F^ = m^F^+n‘ 


= n^F^+m^F^ 
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= -2mn(FpF2) 

= m F^^+n^n (fgg+2F^2)+'' ^22 


= '’^^44 

= 4mV(F^^+F22-2F^2^*(>"^-"^^%6 

= m[2m^{F^^-F^^)-{m^-nh{2Fy^-F^^)] 

= -mn[2n^(F^pF22)+(i"^-n^)(2F^2~f^5g)] 
= -2mn(F^j-F23) 


= mn(. 44-F55) 

= mV(F^^>F,2-Fgg)4^(m^+n^)F^2 

= m’F^3+n^F23 
n F^3+m F^j 
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APPENDIX D 
MATERIAL PROPERTIES 

This appendix contains all the material values for Thomel 3'H) 
Graph ite/Narmco 5208 used in this study. All data were presented in 
reference [22] except for the 0 values in Table 0.2, which came from 
reference [23]. Fig. D.l represents the stress-strain response for 
T300/5208. Table D.l contains the Ranri>erg-0sgood coefficients and 
strength parameters used while Table 0.2 contains the thermal/ 
moisture properties of T300/5208. 

In Table D.l, a* corresponds to the stress at which the Ramberg- 
Osgood coefficients and become applicable. The values of the 
Tsai-Wu interaction terms, ^23* equal to 

-0.58xl0~^^/(PSI)^, the value given in reference [18]. No data was 
found showing the hygrothermal variation of the interaction terms. 
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Fiqure D.l. Stress-Strain Curves for T.100/S?08 Graphite Fnoxy 



TABLE 0.1 Ramberg-Osgood Parameters for Graphite- Epoxy T300/5208 


Curve 

Elastic 

Modulus 

(MSI) 

Elastic 

Limit 

(KSI) 

"l 

-n, 

K^(PSI ') 

0^ 

(KSI) 

"2 

K^(PSI ^) 

m 

X 

X 

27.6 

5.0 

2.555 

1.8074x10"^® 

- 

- 

- 


18.9 

mm 

2.357 

3.0861x10"^^ 

- 

- 

- 


1.6 

0.1 

1.728 

1. 8405x10"^’ 

- 

- 

- 


1.8 

1.5 

3.487 


- 

- 

- 


1.6 

0.1 

1.728 

1.8405x10'^^ 

- 

- 

- 


1.8 

1.5 

3.487 

3.9753xl0‘^® 

- 

- 

- 

E9 


0.05 

1.361 

4.2901x10’^^ 

2.82 

4.668 

1.6842x10"^^ 


1.05 

0.05 

1.361 


2.82 

4.668 



1.05 

0.05 

1.361 

4.2901x10"^° 

2.82 

4.668 

1.6842x10'^^ 

















































TABLE D.2 Hygrothemal Properties for Graphite -Epoxy T300/5208 


r' 


Property 


t 

XX 

c 

XX 

t 

^'L 

c 

n 

t 

Z2 

C 

zz 

'n 

^2 

xy 


Percent Retention of Room Temp., OK Moisture Property 


iRoom Temp, j • 

iOX Moisture Temperature ••70* F I Temperature -260* F i Temperature -350* F 


I -Elastic 
; Modulus- 
, (Msi) 


-t- 


i - 


27.6 

18.9 

1.6 


Percent Weight Gain 


0.00 
100. C 


Percent Weight Gain 


Percent Weight Gain 


100.0 


f 


I 


1.8 
1 .6 
1.8 

1.05 

1.05 
1.05 


100.0 I 88.2 

- ♦ - --- 

100.0 ; 

'' 100.0 ' 88.2 

! - t 

ilOO.O i 




■ 100.0 

100.0 


1100.0 




TABLE 0,2 continued 




TABLE D.2 continued 
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Room Temp. 
jo% Moisture 

1 -V- 

i 

1 

Percent Retention of Roori Temp. , 0% Moisture Property 

Temperature -70® F 

Temperature -260® F 

Temperature -350^ F 

Percen 

0.00 

J 

t Weight 
0.83 

L 

Gain 

1.13 

Percen 
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t Weight 
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: Gain 
1.13 
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Gain 

1.13 

i 0.38 
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. - 4_.__ . H 

89.5 
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- 

- 
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- 

- 
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- 
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- 

- 
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- 
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1 
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- 
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TABLE 0.2 continued 


Room Temp. 

i 0^ Moisture 

Property I -Temperature 

1 Coefficient- 

' (uin/1n/°F) 

? 

( 

Temperature Coefficient (pin/in/°F) 

Temperature -70® F 

Temperature -260° F 

Temperature -350° F 

L _ _ 

Percen 

0.00 

t Weighi 
0.83 

Gain 

1.13 

Percen 

0.00 
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t Weight 
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1,13 

Ti 1 0.01 

X ! 

0.01 

- ...J 

- 

- 

0.14 

- 

- 

0.40 

- 


i 11,0 


- 

- 

17.8 

* 

, 1 


21.0 

- 


CM 

O 

11.0 

'_J 

- 

17.8 

[lu 

- 

21.0 

- 

- 



TABLE D.2 continued 



Room Temp. 
0 % Moisture 

• Ma 4 C f 1 1 1TA 

Moisture Coefficient {in/1n/*Wt Gain) 

Prnnpr*t\y 

Temperature -70* F 

L, 

Temperature -260° F 

Temperature -350° F 

• • P vy p 1 ^ VU 1 ^ 

1 Coefficient- 
I (in/1n/m) 

Percer 

0.00 

^t Weigl 
0.83 

Gain 
1 13 

Percen 

0.00 

t Weigh! 
0.83 

: Gain 
1.13 

Percen 

0.00 

t Weight 
0.83 

: Gain 
1.13 

B 

h- " -i 

0.0 

^0 

0.0 

0.0 

- 

- 

- 

- 

- 

- 


0.0 

0.0 

0.0049 

0.0061 


- 

- 

- 

- 

- 

- 

0.0 

0.0 

0.0049 

1 

0.0061 

- 

- 

j 

- 

- 

- 


no 

r>o 
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APPENDIX E 

INPUT MODIFICATIONS FCfl NONCOM2 


The modifications to NONCOMl reqL'ire two changes to the users 
guide given in reference [3]. One is the input of the Tsai-Wu inter- 
action terms F^^, F^^ and with the uniaxial strengths. The 

other change is the simultaneous input of the thermal and hygroscopic 
properties rather than separately as before. 


E.l The following replaces card 27 of reference [3]. 


Card 27 (6E12.6) 
Column 

1-12 

SL23(K) 

U-24 

SL13(K) 

25-36 

SL12{K) 

37-48 

XF23(K) 

49-60 

XF13(K) 

61-72 

XF12(K) 


Contents 

= Ultimate stress for T 23 - Y 23 
= Ultimate stress for - Y'J 3 
= Ultimate stress for t.j 2 - y -|2 
* F 23 interaction term 
= F ^2 interaction term 
= F ^2 interaction term 


E.2 The following cards replace cards 29 to 50 of reference [3] and, 
when hygrothermal analysis is required, are ■’‘nserted between cards 
28 and 51 as given in reference [3]. 

The following cards are repeated NDIFM time (for each different 
material ). 

K=l, NDIFM 
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Card 29 (1615) 

Col utnn Contents 

1-5 NTE11(K,1) = Number of linear segmented temperature 

points for tensile modulus , ercent 
retention curve 

6-l<^ NME11(K,I,1) = Aumber of linear segronted moisture 

points for F-j^ tensile modulus percent 
retention curve at Ith temperature 

etc. I=1,NTE11(K,1) 

The following two cards are repeated N'Zn(K,l) time (1=1, 
NTEll(k,D). 

Card 30 (1E10.3) 

Column Cr ntonts 

1-10 TMPE11(K,I,1) = Tempetat'^re at Itii temperature 

Car^. 31 (8E10.3) 

Col umn Contents 

1-10 PME11(K,I,J,1) = Moisture content at Ith temperature 

11-20 PRDE11(K,I,J,1) - Percent retention of E^^ tensile 

modulus at Ith temperature and Jth 
moisture content 


ere. 


repeated J=1 ,NME11(K,I,1 ) 



Contents 


Card 32 

(1615) 

Column 

1-5 

NTE11{K,2) 

6-10 

NME11(KI,2) 

Card 33 

(1E10.3) 

Column 

1-10 

TMPE11(K,I,2) 

Card 34 

(8E10.3) 

Column 

I-IO 

PME11(K.I,J,2) 

11-20 

PRDl*11(K,I,J,2) 

Card 35 

(16T5) 

Column 

1-5 

NTE^2(^,1) 

6-10 

NME22(K,I,1) 

Card 36 

(1E10.3) 

Col umn 

1-10 

TMPE22(K,I,1) 

Card 37 
Column 

1-10 

PME22(K,I,J,1) 

n-20 

PRDE22(K,I.J,1) 


= I Same as cards 29 to 31 but for ^ 
compressive modulus 

Contents 


1 

*| Same as cards 29 to 3l but for E 22 
tensile modulus 
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Card 38 (1615) 

Col man 

1-5 NTE22(K.2) 

6-10 NHE22(K.I,2) 

Card 39 (1E10.3) 

Column 

1-10 WE22(K,I,2) 

Card 40 (8E10.3) 

Column 

1-10 PME22(K,I,.^2) 

11-20 PR0E22(K.I,J.2) 

Card 41 (1615) 

Column 

1-5 NTE33(K,1) 

6-10 NME33(K.M) 

Card 42 (1E10.3) 

Column 

1-10 TMPE33(K.I.l) 

Card 43 (8E10.3) 

Column 

1-10 PMC33{K.I,J.l) 

11-20 PR0E33(K.I,J.l) 


Contents 

-I 


A 

) 


Same as cards 29 to 31 but for E ^2 
compressive modulus 




Contents 


Same as cards 29 to 31 but for E^^ 
tensile modulus 


=1 
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Card 44 {1615) 



Coluiw 


C^n^nts 


1-S 

NTC33(K,2) 

\ 


6-10 

NHF331K.l,2) 

\ 


Card 4> 

(ir.10.3) 



^q^lumn 


\ 


1-10 

T>irE33(K.l.2) 

[ Same as cards 29 to 31 

but for E 



1 


Caj^ 46 

(8fl0.3) 

compressive modulus 


Colwnn 




1-10 

PME33(K.I.J.2) 

-1 


11-20 

PR0E33(K.1»J.2) 

= 1 


CvW'd 47 

{16151 



Column 


Contents 


1-S 

NTG23(K) 

\ 


6- 10 

NMC,23(k.l) 

1 


Card 48 

(lflO.3) 



Column 


\ 


1-10 

TMf’r,23(K.I) 

’ 1 Same as cards 29 to 31 

but for 

2j 

Card 48 (SnO.3) 

shear modulus 


Column 




1-10 

PMC.23(K,l.vl) 

1 


11-20 

prog:.?! Ik. 1. 0) 

1 




Contents 


\ 2 <* 


Card 50 

(1615) 

Column 


1-5 

NTG13(K) 

6-10 

NMG12(K.I) 

^rd 51 

(lt\0.3) 

(^Iw 


1-10 

■mP613(K,l) 


•E10.3) 

Co Ur,. 


1-lC 

PMG13(K,I.J) 

11-20 

PRD613(K.I,J) 

Card 53 

(1615) 

Column 


1-5 

NTG12(K) 

6-10 

f#«12(K.I) 

Card M 

(1E10.3) 

Column 


1-10 

THPG12(K,I) 

Ca rd 55 

(8E10.3) 

Column 


1-10 

PHG12(K.I.J) 

11-20 

PRDGi:(K.I,J) 



- \ SaM as cards 29 to 31 but for 

\ 

sheir no«^!i 


-1 

•1 


Contents 

■I 

•1 


Same as cards 29 to 31 but for 
shear modulus 


■I 

■1 
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Cajr^ 

(16151 



Cqlim 


\ 

Contents 

— *■ 

1-5 

NTU23(K.l) 

“ i 

1 

6-10 

NMU23(K.1,11 

. 1 

1 

Card $7 

(1E10.3) 



Column 




\-\0 


\ 

Sane as cards 29 to 31 but for 

Card 58 

(8E10.3) 

tensile Poisson's Ratio 

Column 




1-10 

PMU23(K, 1.0,11 

* i 

1 

11-20 

PR0U23(K, 1.0.1) 

1 

X 


Card 59 

(1615) 



Column 



Contents 

1-5 

NTU23(S ,2) 

-I 


6-10 

NHU23(K.I.2) 


1 

Card 68 

(1E10.3) 



Cojutnn 




1-10 

TMPU23(K.I.2) 

. 1 
1 

1 

‘ Same as cards 29 to 31 but for 

I 

Card 61 

(8E10.3) 

compressive Poisson's Ratio 

Column 




I-IO 

PMR23(K. I. 0.2) 

■1 


11-20 

PROU23(K, I. 0,2) 





Card 62 

(1615) 


Column 


Contents 

V 

1-5 

NTUU(K.l) 

= 1 

6-10 

I«U13{K,I,1) 


Card 63 

(1E10.3) 


Column 


\ 

1-10 

TM>U13(K.U1) 

= 1 Sane as cards 29 to 31 but for 
1 

Card 64 

(8E10.3) 

tensile Poisson's Ratio 

Column 



1-10 

PMU13(K.I,J,1) 

\ 

11-20 

PRUU13(K,1,J.1) 


Card 65 

(1615) 


Column 


Contents 




1-5 

NTU13(K,2) 

■ 1 

6-10 

NMC13(K,U2) 

M 

Card 66 

(1E10.3) 


Column 


\ 

1-10 

TMPU13(K,I.2) 

= 1 Same as cards 29 to 31 but for 

C aj d 6 7 

(8E1Q.3) 

compressive Poisson's Ratio 

Column 



1-10 

PMU13(K.I.J,2) 


11-20 

PRDU13(K.I.J,2) 

- 1 



Card 68 (1615) 
Column 


Contents 


1 o2 


1-5 

NTU12(K.l) 

•1 

6-10 

NHU12(K.I.1) 

= ) 

Card 69 (1E10.3) 


Column 


\ 

1-10 

TMPU12(K.I.l) 

= J Same as cards 29 to 31 but for 
tensile Poisson's Ratio 

Card 70 (8E10.3) 

Column 



1-10 

PMU12(K.I.J,1) 


11-20 

PR0U12(K,1,J.1) 

= 1 

Card 71 

(1615) 


Column 


Contents 

\ 

1-5 

NTU12(K,2) 

= 1 

6-10 

NMU12(K.I,2) 

= ) 

Card 72 (1E10.3) 


Column 


\ 

1-10 

TMPU12(K,I,2) 

= 1 Same as cards 29 to 31 but for 
compressive Poisson's Ratio 

Card 73 (8E10.3) 

Column 



1-10 

PMU12(K,I.J,2) 

= 1 

11-20 

PR0U12(K,I.J.2) 

■ 1 
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Card 74 (1615) 

Column 

1-5 NTS^i(K.l) 

6-10 NMSll(K.I.l) 

Card 75 (1E10.3) 

Column 

1-10 TMPSll(K.I.l) 

Card 76 (8E10.3) 

Column 

1-10 PMSlKK.I.J.l) 

11-20 PRDSll(K.I.J.l) 

Card 77 (1615) 

Column 

1-5 NTS11(K,2) 

6-10 NMS11(K.I.2) 

Card 78 (1E10.3) 

Column 

1-10 TMPS11(K,I.2) 

Card 79 (8E10.3) 

Column 

1-10 PMS1KKJ.J.2) 

11-20 PRDS11(K,I.J,2) 



Contents 


= j Same as cards 29 to 31 but for 
tensile strength 

Contents 


* Same as cards 29 to 31 but for X 

I ' 

compressive strength 
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Card 80 (1615) 


Column 


Contents 

1-5 

NTS22(K,1) 

= 1 

6-10 

NMS22(K,I,1) 


Card 81 

(1E10.3) 


Column 


\ 

1-10 

TMPS22(K,I.l) 

= ^ Same as cards 29 to 31 but for 
1 

Card 82 

(8E10.3) 

tensile strength 

Column 



1-10 

PMS22(K.I.J.l) 


11-20 

PRDS22(K.I.J.l) 

■) 

Card 83 (1615) 


Column 


Contents 

1-5 

NTS22(K,2) 

-1 

6-10 

NMS22(K,I,2) 

•) 

Card 84 (1E10.3) 


Column 


\ 

1-10 

TMPS22(K,I,2) 

= 1 Same as cards 29 to 31 but for 

Card 85 (8E10.3) 

compressive strength 

Column 


\ 

1-10 

PMS22(K,I.J.2) 

=1 

11-20 

PRDS22(K.I,J,2) 
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Card 86 

(1615) 


Column 



1-5 

NTS33{KJ) 

■1 

6-10 

N^•.S33(K.U1) 


Card 87 

(1E10.3) 


Column 



1-10 

TMPS33(K»M) 

- \ 
1 

Card 88 

(8E10.3) 

1-10 

PMS33(K.I»J.l) 

= 1 

11-20 

PR0S33(KJ.J.l) 


Card 89 

(1615) 


Column 



1-5 

NTS33(K,2) 

1 

6-10 

NMS33(K.U2) 


Card 90 

(1E10.3) 


Column 



1-10 

TMPS3o(K.I,2) 

- ( 
\ 

Card 91 

(8E10.3) 

\ 

1-10 

PMS33(K.I.J,2) 

= 1 

11-20 

PRDS33(K,I.J,2) 



Contents 


Same as cards 29 to 31 but for 
tensile strength 


Contents 


Same as cards 29 to 31 but for 
compressive strength 
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Card 92 (1615) 

Column 

1-5 NTS23(K) 

6-10 NMS23(K.I) 

Card 93 (1E10.3) 

Column 

1-10 TMPS23(K.I) 

Card 94 (8E10.3) 

Column 

1-10 PMS23(K,I.J) 

11-20 PRDS23(K,I.J) 

Card 95 (1615) 

Column 

1-5 NTS13(K) 

6-10 NMS13(K.I) 

Card 96 1E10.3) 

Column 

1-10 TMPS13(K,I) 

Card 97 (8E10.3) 

Column 

1-10 PMS13(K.I,J) 

11-20 PRDS13(K,I,0) 


Contents 

-I 


I 


Same as cards 29 to 31 but for $23 
shear strength 


•I 
= ) 



Contents 


= ; Same as cards 29 to 31 but for S,^ 

I 

shear strength 

■I 
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Card 98 (1615) 


C''lumn 


Co.itents 

1-5 

NTS12(K) 


6-10 

NMS12(K,I) 

■) 

Card 9S 

(IEIO.3) 


Column 


\ 

1-10 

TMPS12(K,I) 

= > Same as cards 29 to 31 but for S^2 

Card 100 (8E10.3) 

/ 

shear strength 

Col uinn 



1-10 

PMS12(K,I,J) 

-1 

11-20 

PRDS12(K,I,J) 


Card 101 

(1615) 


Column 


Contents 

1-5 

NTF23(K) 

-1 

6-10 

NMF23(K,I) 


Card 102 

(1E10.3) 


Column 


\ 

1-10 

TMPF23(K,I) 

= 1 Same as cards 29 to 31 but for F23 

Card 103 

(8E10.3) 

interaction term 

Column 



1-10 

PMF23(K,I,J) 

= 1 

11-20 

PRDF23(K,I,J) 

= 1 
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Card 104 (1515) 

Column 

1-5 NTF13(K) 

6-10 NMF13(K,I) 

Card 105 (1E10.3) 

Col umn 

1-10 TMPE13(K,I) 

Card 106 (8E10.3) 

Column 

1-10 PMF13(K.I.J) 

11-2C PRDF13(K,I,J) 

Card 107 (1615) 

Column 

1-5 NTF12(K) 

6-10 NMF12(K,I) 

Card 108 (1E10.3) 

Column 

1-10 TMPF12(K,I) 

Card 109 (8E10.3) 

Column 

1-10 PMF12(K,I.J) 

11-20 PRDF12(K,I,J) 



Contents 


= } Same as cards 29 to 31 but for F,- 
I ij 

interaction term 

-) 

Contents 



Same as cards 29 to 31 but for F^^ 
interaction terms 
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Card no (1615) 



Column 



Contents 

1-5 

NTALl(K) 

= 1 


6-10 

NMALI(K.I) 

f 

Same as cards 29 and 30 but for 

Card 111 

(1E10.3) 

temperature coefficient 

Co 1 umn 




1-10 

TMPAL1(K,I) 

A 

) 


Card 112 

(8E10.3) 


Column 



Contents 

1-lC 

PMAL'i(KJ.O) 


Moisture content at Ith tempearture 

11-20 

pri'Ali(k:.i,j) 


temperature coefficient at Ith 
temperature and Jth moisture content 

etc. 

repeated J=1 , 

NMALl(K.I) 


Card 113 (1615) 



Column 



Contents 

1-5 

NTAL2(K) 

= 1 


6-10 

NMAL2(K,I) 

= 1 


Card 114 

(1E10.3) 



Column 


t 


1-10 

TMPAL2(K,r) 

Same aS cards 110 to 112 but for 

Card 115 

(8E10.3) 

temperature coeffic’ent 

Column 




1-10 

PMAL2(K,I,J) 

= 1 


11-20 

PR0AL2(K,I,J) 
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Card 116 (1615) 

Column 

1-5 NTA5.3;.v) 

6-10 NMAL3(K,I) 

Card 117 (1E10.3) 

Column 

1-10 TMPAL3(K,I) 

Card 118 (8E10.3) 

1-10 PMAL3(K,I,J) 

11-20 PRDAL3(K.I.J) 

Card 119 (1615) 

Column 

1-5 NTBTl(K) 

6-10 NMBT1(K,I) 

Card 120 (1E10.3) 

Column 

1-10 TMPBT1(K,I) 

Card 121 (8E10.3) 

1-10 PMBTl(K.I.J) 

11-20 PRDBT1(K,I,J) 


Contents 

A 

-\ 


= I Same as cards 110 to 112 but for 
temperature coefficient 

A 

A 

Contents 


i 


A 

=) 


Same as cards 110 to 112 but for 
moisture coefficient 



Card, \:2 (IbIS) 

Column 

1-5 ntbt:’(k' 

6-10 nmbt:{k,i) 

Card 123 (IE 10. 3) 
Column 

1-10 THPBT2(K.I) 

Card 124 (8fl0.3) 

Cq^lumr 

1-10 PMBT2(K.l,J) 

11-20 PR08T2(K.I.J) 

Card 125 (1615) 

Colunm 

1-5 NTBT3(K) 

6-10 NMBT3(K.I) 

Card 126 (U 10.31 
Coluirai 

1-10 1MPBT3(K.I) 

Cai'd 127 (Bn0..11 
Colunm 

l-U) PMBT3(k.l.Jl 

11-20 PROBT.UK.l.O) 


Center. .•; 

I 

) 


Same as cards 110 to 112 bat for 
moisture coefficient 


\ 

\ 

Contents 

= I 
) 


Same as cards 110 to 112 bat for 
moisture coefficient 


I 

I 



